Nonlinear dynamics of polariton scattering in semiconductor microcavity: bistability 

vs stimulated scattering 
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We demonstrate experimentally an unusual behavior of the parametric polariton scattering in 
semiconductor microcavity under a strong cw resonant excitation. The maximum of the scattered 
signal above the threshold of stimulated parametric scattering does not shift along the microcavity 
lower polariton branch with the change of pump detuning or angle of incidence but is stuck around 
the normal direction. We show theoretically that such a behavior can be modelled numerically by 
a system of Maxwell and nonlinear Schrodinger equations for cavity polaritons and explained via 
the competition between the bistability of a driven nonlinear MC polariton and the instabilities of 
parametric polariton-polariton scattering. 

PACS numbers: 71.36.+C, 42.50.p, 42.65.k, 42.65.Sf, 68.65.+g 



The effect of parametric scattering of semiconductor 
microcavity (MC) polaritons has been established re- 
cently 0, fi 1 H 1 i i i [3 [H [13 and at- 
tracted a significant interest. Extremely large parametric 
gains ilfl] and striking transformations of the polariton 
parametric luminescence with the increase of intensity 
of resonant coherent excitation |^ |M ISB- have been 
demonstrated and analyzed theoretically 0,0 IE 0| ■ 
However, only the case of specific pump conditions known 
as magic angle excitation, when the parametric scatter- 
ing threshold intensity is minimum, has been investigated 
so far in detail in case of stationary excitation. 

In this Letter we investigate the parametric MC 
polariton-polariton scattering varying the cw pump an- 
gle and energy away from the magic angle conditions. 
According to the existent theory 0,M 13 > the maxima of 
scattered signal have to shift along the blueshifted and 
renormalized lower polariton (LP) branch. To the con- 
trary, we found experimentally that above the threshold 
of stimulated parametric scattering the signal maximum 
stays approximately in the normal direction, notwith- 
standing the change of the excitation conditions 0| . 

We show numerically that a qualitatively similar be- 
haviour can be demonstrated by a system of Maxwell 
and nonlinear Schrodinger equations for MC quantum 
well excitons, accounting for the exciton-cavity photon 
coupling into MC polariton and polariton-polariton scat- 
tering. In the real life, many other scattering pro- 
cesses can play a role, e.g., polariton-phonon scatter- 
ing polariton- free carrier scattering But it 
turns out that the polariton-polariton scattering alone 
demonstrates already sharp transformation of the scat- 
tering process, which results in the jump of the scattered 
signal from the oblique to normal direction, in a qualita- 
tive agreement with the experiment. 

In the part of polariton-polariton scattering, our de- 
scription is analogous to that in 0, 0, 13 ■ As in 
we neglect the saturation term, responsible for the Pauli 
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FIG. 1: The experimental geometry (a); measured energy- 
momentum spectra of the LP photoluminescence under a res- 
onant excitation at hflo = 1452.7 meV and i9 = 12.5° with 
intensity above the scattering threshold (b). 

principle in exciton-photon interaction. It can be shown 
that its influence is not strong in case of MC with our 
parameters. However, one difference is important. We 
treat the amplitude of the driven polariton mode not as 
a fixed external parameter, but as an internal parameter 
established in the nonlinear system under a fixed exter- 
nal pump. Well known in the theory of driven nonlinear 
optical cavities without a polariton resonance (see, e.g., 
in ^3)' this description allows to explain bistabil- 
ity, symmetry breaking, pattern formation, bifurcations 
to chaos. As we show in this Letter, the bistability can 
play an important role in the nonlinear kinetics of MC 
polariton scattering as well. 

Experimentally, we have investigated photolumines- 
cence (PL) from GaAs/AlAs MC containing six 10-nm- 
thick InGaAs quantum wells (QW's) in 3/2 A GaAs cavity 
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FIG. 2: Lower polariton branches; measured (circles) and 
calculated for |Pol^ = (a), 0.14 (b), and 0.44 (c). Inserts in 
panels b, c show the corresponding dampings. The lines are: 
signal (solid) and idler (dashed). Dash-dot-dot, dash-dot and 
dash lines in panel a are idlers calculated for IVof = and 
different pump angles and detunings (i9, A) = (12°, 0) - magic 
excitation, (14° ,0) and (14°, 0.3 meV). [A = Mo - -ELp(ko).] 
Square, circle, and triangle are the corresponding pump (P) 
locations. Thin lines in panels b, c are signal and idler at 
I Pop = from panel a. 

with cavity mode (cavity mode - QW exciton detuning 
-0.5 meV), see the experimental geometry in Fig. A 
sample was mounted in a cryostat with a wide angular 
access (70°) at a controlled temperature of T=1.8-20 K. 
The experiment was done under excitation with a tun- 
able cw Ti-Sapphire laser under a changeable angle of 
incidence ??o- The angle resolved measurements of the 
MC emission and the use of cw excitation have allowed 
to fix both the energy and wavevector of the emission 
with a high precision. The details on the sample and 
experimental set-up can be found in "4] . 

Figure ^ shows the distribution of the scattered 
light vs. energy fko and in-plane momentum k = 
hcusind/c. The excitation is at do = 12.5° and HHq — 
1452.7 meV (fco ~ 1.6fJ,m.~^), with high excitation power 
of 880 W/cm^, which is above the threshold of stimulated 
parametric scattering. Figure shows two pronounced 
emission maxima at A; ^ and k ^ 2ko. A closer look 
shows that this cannot be explained within the existent 
theory of the MC parametric scattering 0j H 0] j even if 
the LP spectra renormalization is taken into account. 



The bare LP dispersion i?Lp(k) has been determined 
experimentally from the angle resolved PL measurements 
under a weak above band gap excitation, see in Fig. ^p.. 
The inflexion point of the LP branch is at « 12°, which 
corresponds to momentum « 1.5 • 10** cm~^ and energy 
1452.2 meV. 

At low excitation intensity the LP polariton-polariton 
parametric scattering should take place with the en- 
ergy and in-plane momentum conservation, k + k' = 
2ko, i?Lp(k) -I- £'Lp(k') = 2?ii7o. Two sharp emission 
peaks, signal (S) at k^ and idler (I) at k^, are predicted 
at the intersection of the signal and idler branches i?LP (k) 
and 2hflo — i?Lp(2ko — k), see in Figds. Under a weak 
pump into the inflexion point (magic angle) the signal 
peak is predicted at kg ~ 0. With the shift of the pump 
angle and energy away from the inflection point S and I 
peaks shift along the LP branch towards kg < 0, ki > 2ko 
as shown by arrows in Fig. ^p,. 

With increasing the pump intensity the signal and idler 
branches become renormalised via the parametric cou- 
pling 0,01, see in Fig.EJ*. Simultaneosly the dampings of 
signal and idler branches at the peak positions approach 
zero (see insert in panel b). At some critical pump inten- 
sity the damping changes sign (becomes gain) , which sig- 
nals the instability and means the threshold of the para- 
metric stimulated scattering. However, the prediction 
within this theoretical model is still kg < 0, ki > 2kQ, 
whereas on the experiment we observe the signal and 
idler maxima at fcs ~ 0, ki ^2kQ. 

In order to model the polariton scattering dynamics, 
we solve numerically the system of coupled equations for 
£qw, the electric field on the QW, and V{k,t), the exci- 
ton polarisation integrated over the QW thickness, 

[ihf^ - Emc] £Qw{k, t) = a£,^t{k, t) + l3V{k, i), (1) 

[ift^ - Ex]V{k,t)^A£ciwik,t) + at) (2) 

+F Viq, t)V{q\ t)V*{q + q' - k, t). 
q,q' 

Here £cxt — £cxt{t)e~^^^°*'6k,ko is the external electric field 
of the incoming pump far from the MC which is treated 
as coherent macro-occupied mode with fixed H^Iq and 
ko- £cxt{t) is the pump amplitude, changing slowly with 
time. i?Mc(fc) a-iid Ex are empty microcavity and QW 
exciton frequencies. F and A are the QW exciton in- 
teraction constant and susceptibility, a{k) and P{k) are 
the MC response coefficients. In these notations, the MC 
polariton Rabi frequency is y/Aj3. Both Eucik) and Ex 
are complex, with imaginary parts accounting for finite 
lifetimes of the MC mode and QW exciton. ^(fc, t) is 
the stochastic Langevin force. Eucik), a{k), and (3{k) 
are calculated via the scattering matrices of the upper 
and lower MC Bragg mirrors. Ex and A are fitted pa- 
rameters. The energy is measured in meV, the units of 
polarisation and electric field are such that F = 1. The 
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FIG. 3: The calculated time-integrated energy-momentum 
spectra of the scattered light from MC when the pump pulse 
amplitude is slightly below (a) and above (6) the threshold. 
Dotted lines are the signal and idler branches from Fig.|5J) cal- 
culated at the parametric scattering threshold \Vof ~ 0.14. 

in-plane momentum is two-dimensional, but we solve nu- 
merically only a simplified one-dimensional model. 

Equation ^ is the Maxwell equation with the ex- 
citon polarization in the resonance approximation for 
the cavity. Equation ^ is the inhomogeneous nonlin- 
ear Schrodinger equation with damping and two types of 
sources: coherent external pump and stochastic Langevin 
noise. Eas. Hll2ll allow to describe the influence of coher- 
ent pump and quantum fluctuations simultaneously in a 
very convenient for numerical modelling way. 

The numerical solution of Eqs. (|1I2|) for long 1 ns) 
nearly rectangular excitation pulses demonstrates a sharp 
threshold-like transition from a picture agreeing with the 
described above stationary model to a completely differ- 
ent one, compare Figs-OJa and b. Note that the difference 
in the pump intensities between panels a and & is « 1 % 
only. The averaged intensity of scattered noise increases 
above transition by several orders of magnitude (note the 
logarithmic vertical scale in Fig.|3|). 

Figure 01 explains the time kinetics of the scattering 
above the threshold pump intensity. The shape of the 
excitation pulse f g^t (dash-dotted line in Fig. 2Ji) was 
chosen in such a way that the threshold intensity is ap- 
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FIG. 4: The time dependences: of input excitation pulse, 
calculated exciton polarisation and absorption coefficient (a); 
of angular spectra of the scattered light (6). Exciton polarisa- 
tion vs pump intensity during the excitation pulse (c). Insert 
shows the magnified part of the region of transitions. The 
time delay in panel c between red, black, and triangled dots 
is 1, 10, and 100 ps, respectively; the triangles are labelled by 
time in ps. Solid S-shaped curve is the solution of Eq. Q. 



proached slowly during the pulse duration. Actually two 
sharp transitions take place (at t « 600 and 700 ps). 
Both transitions are characterized by the jumps of the 
driven mode amplitude Vq itself, see solid line in Fig.^. 

In order to understand these transitions, we have to 
investigate the stability of solutions of Eqs. H1I2|I in 
the case of stationary external fleld £e-xt{t) ~ const 
with only one macroscopically flUed mode, i.e., of the 
form V{k,t) = V{k,t) + 4,fcoPoe-*^°*, and £Qw{k,t) = 
£{k,t) + Sk,ko£oe-'^°*. Here {V/Vol \£/£o\ < 1- In zero 
order we get the following cubic equation for the ampli- 
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tude of the driven mode Vq, The influence of the bistability was not analyzed in the 

existent theories of the MC polariton scattering 0, 0, S 

[Ao{Ao + Axc)-Af3]Va-AoF\VafVo^Aa£,^t, (3) 
where Aq = h^o - EMc{ko),^xc = Euciko) - Ex- 

Solution of Eq. Q is shown as the sohd S-shaped hne As to the stability conditions, Hnearising Eqs. (UEl, 

in Fig. H:. The S-shape is well known in the theory of we have to solve the linear eigenproblem for parametri- 

nonlinear cavities (see, e.g., in Refs. [3113). brings cally coupled signal £^V{k,t) = £,P{k)e-"^^ and idler 

a bistability into the behaviour of the MC polaritons. £*, :P*(fc, t)e-2*"''* = £,r{k)e-''^* (where k ^ 2ko - fc). 
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and check the sign of the imaginary parts of the eigenen- 
ergies. Equation |@J appears to be analogous to that 
discussed in detail in Refs. 0,1101. We use the MC-X 
basis instead of the UP-LP one. There is an advantage 
of the MC-X basis, because the dependence of the inter- 
action constant F on k can be neglected. Then, in the 
space representation, it corresponds to a local contact 
interaction F J \'P{x)\'^dx, which allows to write a very 
efficient numerical code for Eqs. (|1I2I) . 

The idler and signal branches discussed above [see in 
Figs.|3)|3] are the eigenvalues of Eq. Q calculated with 
Vq — 0.14, a threshold value for the stimulated scattering 
into ks < 0, ki > 2kQ. And indeed. Fig. shows that 
at i < 600 ps the scattered noise is maximum at the pre- 
dicted angles kg < 0, k^ > 2ko. However, at t « 600 ps, 
instead of the stimulated scattering into these modes, an- 
other instability develops. The trajectory of the system 
on the [[^'(/co, |fGxt(t)P] plane, see in Fig.^, shows 
that the instability at t « 600 ps is the jump of the ko 
mode between the lower and upper branches of the S- 
shaped curve. 

In the empty cavity with quadratic dispersion, the up- 
per S' branch is usually stable 0, ^ with 
inflective LP dispersion, it can become unstable against 
parametric scattering. It can be seen from, e.g.. Fig. 12:;, 
showing the signal and idler branches calculated on the 
upper part of S-curve at jPoP = 0.44. Instead of damp- 
ing, large gain is realized for polariton modes just with 
ks ~ 0, ki ^ 2ko. As a result, a stimulated scattering 
into these modes develops (see in Fig. 01*). It is feed up 
by the increased absorption, see dashed line in Fig. ^Ji. 
But the pumped mode becomes unstable and eventually 
jumps back into the lower position (at t « 700 ps). Now, 
because we already have a well developed scattered noise 
in the system, this lower position is more or less sta- 
ble, although subject to noisy fluctuations because of the 



already developed scattered states and the competition 
between them and the pumped mode. 

We would like to note that within our theoretical model 
the instability of the upper S' branch takes place in a 
narrower region of parameters than in the experiment. 
E.g., for the pump at 1)0 — 14° and 12.5° (as in the 
experiment in Fig. QI)) it occurs for detuning 0.1 < A < 
0.6 meV and 0.15 < A < 0.3 meV, respectively. However, 
the other scattering mechanisms, e.g., with phonons and 
free carriers, which are neglected in our approach, can 
increase substantially the region of the instability. 

To conclude, we have found experimentally that above 
the MC parametric scattering threshold the signal and 
idler maxima do not shift along the renormalised LP 
branch but jump to kg ^ 0, ki ^ 2ko directions. We 
show that qualitatively this behavior can be explained 
via an interplay between the the bistability of the driven 
mode itself and the stimulated parametric scattering. 
The bistability, although well known in empty nonlinear 
cavities, was not accounted for in the existent theoretical 
models of the MC polaritons parametric scattering. It 
acquires a new feature in MC with polariton resonance: 
instability of the upper branch of the S-shaped depen- 
dence of pumped mode polarisation vs pump. Although 
the exciton-exciton scattering may be not the only mech- 
anism causing the instability, this feature is important 
and has to be taken into account. 
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